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ABSTRACT 



We present i?-band photometry of the X-ray transient and candidate black hole binary XTE 
J1650-500 obtained between 2003 May and August with the 6.5m Clay Telescope. A timing 
, analysis of these data reveals a photometric period of 0.3205 ± 0.0007 days (i.e. 7.63 hr) with 

' a possible alias at 0.3785 days (9.12 hr). Our photometry completely rules out the previously 

published spectroscopic period of 0.212 days (5.09 hr). Consequently, we reanalyzed the 15 
archival ESO/VLT spectra (obtained 2002 June by Sanchez-Fernandez et al.) that were the basis 
of the previously published spectroscopic period. We used a "restframe search" technique that is 
well suited for cases when the signal-to-noise ratio of individual spectra is low. For each of roughly 
1.1 million binary ephemerides, we summed all of the spectra in a trial restframe of the secondary 
star, and each restframe spectrum was cross-correlated against a template spectrum. We then 
searched for the set of orbital parameters that produced the strongest cross-correlation value. 
The results confirmed the photometric period of 0.3205 days, and rule out the alias period near 
0.38 days. The best value for the velocity semiamplitude of the companion star is K2 = 435 ± 30 
km s~^, and the corresponding optical mass function is f{M) = 2.73 ± 0.56 M©. The spectral 
type of the companion star is not well constrained because we only have six template spectra 
available to us. The K4V template provides the best match; next best matches are provided 
by the G5V and K2III templates. We also find that the accretion disk dominates the light in 



^Based on observations made with the Magellan 6.5m Clay telescope at Las Campanas Observatory of the Carnegie Institution 
and on observations made with European Southern Observatory (ESO) Telescopes at the Paranal Observatory under programme 
ID 69.D-0644(A) 
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the i?-band where the disk fraction is 80% or higher, although this value should bo treated with 
caution owing to the poor signal-noise-ratio and the limited number of templates. The amplitude 
of the phased i?-band light curve is 0.2 magnitudes, which gives a lower limit to the inclination of 
50 ± 3° in the limiting case of no contribution to the i?-band light curve from the accretion disk. 
If the mass ratio of XTE J1650-500 is similar to the mass ratios of other black hole binaries like 
A0620-00 or GRS 1124-683 (e.g. Q > 10), then our lower limit to the inclination gives an upper 
limit to the mass of the black hole in XTE J1650-500 of Mi < 7.3 Mq. However, the mass can 
be considerably lower if the i?-band flux is dominated by the accretion disk. For example, if the 
accretion disk does contribute 80% of the flux, as our preliminary results suggest, then the black 
hole mass would be only about 4:Mq. 

Subject headings: binaries: spectroscopic — black hole physics — stars: individual (XTE J1650-500) — 
X-rays: stars 



1. Introduction 

X-ray novae (XN) provide the strongest evi- 
dence for the existence of stellar mass black holes. 
XN are interacting binaries where a neutron star 
or a black hole accretes matter from a compan- 
ion (usually a dwarf-like K or M star). The ac- 
cretion rate on to the black hole increases sub- 
stantially during the outburst phase, and hence 
the X-ray luminosity can vary by large amounts 
(e.g. factors of 10^ or more). The majority of XN 
spend most of their time in a "quiescent" state 
where the X-ray luminosity is on the order of the 
optical luminosity of the companion star. It is 
in quiescence where the companion star can be 
best studied in the optical, where the observed 
radial velocity and light curves of the companion 
lead to dynamical mass measurements for the com- 
pact primary. There are 18 X-ray binaries (15 of 
them XN) where the mass of the primary has been 
shown to exceed the maximum mass of a stable 
neutron star (f« 3Mq, Kalogera & Baym 1996), 
confirming the presence of black holes in these sys- 
tems: GRO J0422-F32 (Orosz & Bailyn 1995, Fil- 
ippenko, Matheson, & Ho 1995); A0620-00 (Mc- 
CUntock & RemiUard 1986); GRS 1009-45 (Fil- 
ippenko et al. 1999); XTE J1118-F480 (McClin- 
tock et al. 2001, Wagner et al. 2001); GS 1124-683 
(RemiUard, McClintock, & Bailyn 1992); 4U 1543- 
47 (Orosz et al. 1998); XTE J1550-564 (Orosz et 
al. 2002); GRO J1655-40 (Bailyn et al. 1995); GX 
339-4 (Hynes et al. 2003); H1705-250 (RemiUard 
et al. 1996); SAX J1819.3-2525 (Orosz et al. 2001); 
XTE .11859+226 (Filippenko & Chornock 2001); 
GRS 1915-1-105 (Greiner, Cuby, & McCaughrean 
2001); GS 2000-1-25 (Casares, Charles, & Marsh 
1995); GS 2023-^338 (Casares, Charles, & Naylor 



1992); Cyg X-1 (Gies & Bolton 1986); LMC X- 

3 (Cowley et al. 1983); LMC X-1 (Hutchings ct 
al. 1987). These sources open up the possibility 
of studying general relativity in the strong field 
regime. For example, the stiidy of high frequency 
quasiperiodic oscillations in the X-ray light curves 
of certain XN may lead to a measurement of black 
hole spin (e.g. McClintock & RemiUard 2004 and 
cited references). We must press hard to obtain 
further observations of black hole masses in order 
to fully pursue these opportunities. 

XTE J1650-500 (hereafter J1650) was discov- 
ered by RXTE on 2001 September 5 (RemUlard 

2001) and subsequently reached a peak X-ray in- 
tensity of 0.5 Crab. Based on subsequent obser- 
vations, J1650 was established as a strong black 
hole candidate based on its X-ray spectrum and 
variability in the X-ray light curve (Markwardt, 
Swank, & Smith 2001; Revnivtsev & Sunyaev 
2001; Wijnands, Miller, & Lewin 2001). The ra- 
dio counterpart was discovered with the Australia 
Telescope Compact Array (ATCA) by Groot et 
al. (2001). Further radio observations sampled 
the behavior of XTE J1650-500 along all its X- 
ray states (Corbel et al. in preparation). We high- 
light two key results obtained during its outburst 
phase. First, RXTE observations during the third 
and fourth weeks of the outburst yielded a strong 
X-ray QPO with an rms amplitude of 5.0 ± 0.4% 
at a frequency of v = 250 ± 5 Hz (Homan et al. 
2003b). Second, XMM-Newton observed a broad, 
skewed emission line due to Fe Ka (Miller et al. 

2002) . Those authors argue that their results im- 
ply the primary is a nearly maximal Kerr black 
hole that is delivering its spin-down energy to the 
accretion flow. 
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The first significant observational program in 
the optical was that of Sanchez-Fernandez et al. 
(2002, hereafter SF2002). They observed J1650 
on the night of 2002 June 10 with the fourth 8.2m 
telescope at the European Southern Observatory, 
Paranal, and reported the following orbital ele- 
ments: an orbital period of P = 0.212 ± 0.001 
days and a velocity semiamplitude of = 309 ±4 
km s~^, resulting in an optical mass function of 
/(M) = 0.64 ± 0.03 Mq. The results in this paper 
contradict these findings. 

In this paper we report the results of our pho- 
tometric study of J1650. A time series analysis 
of our photometry rules out the orbital period re- 
ported by SF2002. Consequently we also report 
herein our rcanalysis of the SF2002 data obtained 
from the ESO archives. We show that the spectro- 
scopic period we derive from these data is consis- 
tent with our photometric period P = 0.3205 days, 
and we go on to determine the orbital elements of 
the system. We outline below our observations 
and reductions, and our analysis techniques. We 
end with a brief discussion of the implications of 
our results regarding the fast QPO observed for 
J1650. 

2. Observations 

2.1. Optical Photometry 

We observed J1650 with the 6.5m Clay tele- 
scope at Las Campanas Observatory of the 
Carnegie Institution 2003 May 31 and June 1 us- 
ing the Magellan Instant Camera (MagIC) and 
an ii-band filter. A total of 47 i?-band images 
with typical exposure times between 500 and 600 
seconds were obtained over the two nights in good 
conditions: the average seeing was ~ 0.6 arcsec 
and it was photometric about 75% of the time. 
M. Holman and P. Schechter kindly provided an 
additional 65 images obtained 2003 August 1 and 
2 using the same instrumentation, and the ob- 
serving conditions were comparable to those just 
described. In order to minimize readout time, 
the Magic camera is read out simultaneously us- 
ing four amplifiers. Thus, bias and flat-fielding is 
done separately for each quadrant of the detec- 
tor. These calibrations and the merging of the 
quadrant images into a single image were per- 
formed using the publicly-available MagIC reduc- 
tion pipeline. Stetson's programs DAOPHOT He, 



ALLSTAR, and DAOMASTER (Stetson 1987, 
1990; Stetson, Davis, & Crabtree 1991; Stetson 
1992a,b) were used to extract the stellar intensi- 
ties and derive the light curve for J 1650 (sec Orosz 
& Wade 1999 for a detailed discussion of the over- 
all procedure used). This suite of codes gives 
robust instrumental magnitudes. DAOPHOT fits 
the point spread function (PSF) for each image 
using several relatively isolated bright stars, and 
ALLSTAR uses the PSF and finds instrumental 
magnitudes for all of the stars on an image si- 
multaneously (local background subtraction is in- 
cluded). Finally DAOMASTER iterativcly solves 
for zero-point shifts in the magnitude scales from 
image to image by essentially using all of the sta- 
ble stars as "comparison" stars. Thus the final 
results are insensitive to changes in the seeing and 
to changes in the sky transparency. Fig. 1 shows 
a finding chart for J1650 made from the MagIC 
data. 

2.2. Optical Spectroscopy 

SF2002 observed J 1650 on the night of 2002 
June 10 with the fourth 8.2m telescope at the 
European Southern Observatory, Paranal, using 
the F0RS2 imaging spectrograph with an MIT 
CCD mosaic, the 1200R grism, the GG435 block- 
ing filter, and a fixed slit width of 0.7 arcsec- 
ond. This instrumental configuration yields a 
wavelength range of 5750-7310 A with a resolu- 
tion of about 0.76 A FWHM. SF2002 obtained a 
total of 15 spectra of J1650 in somewhat marginal 
seeing (typically between 1.0 and 1.5 arcseconds). 
They also obtained the spectra of six bright G- 
and K-type dwarfs using the same instrumental 
configuration. We obtained these data from the 
ESO archive, and used tasks in IRAF^ to perform 
the CCD reductions and to extract the spectra. 
Since the fiat field and wavelength calibration im- 
ages were obtained exclusively during the daytime 
hours at Paranal, it was necessary to correct for 
fiexure in the spectrograph by applying small cor- 
rections (< 0.6 A) to the wavelength scales of the 
extracted spectra using bright night-sky emission 
lines. 



^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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3. Analysis 

The Magellan light curve of J1650 showed con- 
siderable variability, so we phased the data on the 
period reported by SF2002 (P = 0.212 days). To 
our surprise, the light curve phased on that pe- 
riod showed considerable scatter. Consequently, 
we searched the Magellan light curve for period- 
icities using the ELC code with its genetic fitting 
algorithm (Orosz & Hauschildt 2000; Orosz et al. 
2002). About 1,106,000 ellipsoidal models were 
generated and compared with the data using the 
absolute deviation as the merit function: 

where y{xi;a) denotes the model value computed 
at Xi, Hi is the observed quantity at the same xi, 
and (Tj is the uncertainty in j/j. This merit func- 
tion is more robust than the function because 
it is less sensitive to a few outlying points. Mod- 
els were computed using a period range of 0.15 to 
0.50 days. After all of the fits were performed, the 
merit function hypersurface was projected onto 
the trial period axis (see Orosz et al. 2002 for a 
more in-depth discussion of this error estimation 
technique); the resulting periodogram is shown in 
the top of Fig. 2. The absolute deviation has a 
minimum value at a trial period of 0.3205 days, 
and a secondary minimum at a trial period of 
0.3785 days. 

Fig. 3 (top) shows the Magellan light curve 
phased on a period of P = 0.3205 days. There is 
no dip in the absolute deviation near the period re- 
ported by SF2002 (P = 0.212 days), and the light 
curve phased on that period is essentially a scatter 
plot (bottom of Fig. 3). As a check on our results, 
we used the 'pdm' task in IRAF, which is an im- 
plementation of the phase dispersion technique of 
Stellingwerf (1978). A similar periodogram was 
obtained. To estimate the error on our photo- 
metric measurement of the orbital period, we did 
additional fits using ELC's genetic algorithm and 
the normal s-s the merit function. We find a 
1(7 uncertainty of about 0.0007 day, so we adopt 
P = 0.3205 ± 0.0007 day. 

We then analyzed the VLT spectra in an at- 
tempt to understand the disagreement of our pho- 
tometric period with the spectroscopic period re- 
ported by SF2002. We tried to extract radial ve- 



locities using the standard cross-correlation tech- 
nique of Tonry & Davis (1979, implemented in the 
IRAF 'fxcor' task). However, we found that the 
spectra wore very noisy and yielded only a few 
marginal measurements of velocity. We therefore 
used the "restframe" analysis that we developed 
for H1705-250 (Rcmillard et al. 1996). This tech- 
nique is similar to the "skew mapping" technique 
sometimes used for cataclysmic variables (Smith, 
Cameron, & Tucknott 1993; Vande Puttc et al. 
2003), and is quite simple to employ. Suppose 
one has a time series of spectroscopic observations. 
If the spectra arc Doppler shifted to zero veloc- 
ity and coadded using the correct orbital elements 
(e.g. the period P, the semiamplitude K2 and the 
time of maximum velocity To), then the absorption 
features of the companion star will appear at the 
same wavelengths in all of the individual spectra. 
Consequently, the lines in the summed spectrum 
will have a higher signal-to-noise ratio than in the 
individual spectra. If, on the other hand, one uses 
the incorrect orbital elements, then the absorption 
lines from the companion will be at different wave- 
lengths in different individual spectra, and hence 
will be "averaged out" in the resulting summed 
spectrum. 

In our implementation of the restframe anal- 
ysis, we used a FORTRAN program that writes 
IRAF scripts to do the Doppler shifting and the 
spectrum summation. The fxcor task was used to 
provide a measure of how well a restframe spec- 
trum matched the template spectrum. In partic- 
ular, we used the peak cross-correlation value as 
the measure of the goodness-of-fit of the template 
spectrum to the restframe spectrum. The wave- 
length region used for the cross-correlation analy- 
sis was 5907-6261A, 6304-6513A, 66IO-68OOA, and 
7000-7239A. This region excludes telluric lines, 
and the Ha emission line. 

We are interested in three orbital elements: the 
period P, the semiamplitude K2 and the time of 
maximum velocity Pq- This three dimensional pa- 
rameter space was searched by making two dimen- 
sional grids of restframe spectra in the K2 — Tq 
plane at several different values of the period P. 
Specifically, we used periods in the range 0.1500 
to 0.5000 days in steps of 0.0020 days. The rest- 
frame spectra for each trial period were stored in 
a separate subdirectory on disk. For each period, 
we constructed a grid of restframe spectra in the 
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K2 — Tq plane as follows. The range considered 
for K2 was 300-600 km in steps of 2 km s~^. 
For To the range was centered on HJD 2,452,436.6 
and the step size was 0.005 days. The extent of 
the range for Tq was adapted to accommodate the 
trial value of the period. An individual restframe 
spectrum was generated for each set of orbital el- 
ements. In total we generated 1,166,146 restframe 
spectra using this same number of scripts stored 
as separate files. 

Next, we did a cross-correlation analysis on 
the restframe spectra. It was therefore neces- 
sary to find a good template spectrum. To do 
this, we adopted a period of 0.320 days found 
from the photometry and considered the rest- 
frame spectra in the K2 — Tq plane at this period. 
We then cross-correlated these restframe spectra 
against all six of the template spectra in turn and 
determined that the K4V star BS5568 provided 
the best match since it gave the strongest cross- 
correlation peaks. Next, all of the restframe spec- 
tra at all periods were cross-correlated in batch 
mode against the K4V template. We made a pe- 
riodogram by parsing the fxcor log files to get 
the peak cross-correlation value for each restframe 
spectrum. The maximum peak cross-correlation 
value within each period subdirectory was saved, 
resulting in the periodogram shown at the bot- 
tom of Fig. 2. The maximum cross-correlation 
value occurs for a trial period of P = 0.320 days, 
in agreement with the photometric results, which 
thereby establishes the orbital period. Thus, we 
rule out the alias photometric period near 0.38 
days and again rule out the 0.212 day spectro- 
scopic period reported by SF2002 (Fig. 2). We 
also made a periodogram using a template with 
spectral type G5V. Overall, the cross-correlation 
values were much lower and the periodogram was 
much noisier. However, there still was a peak near 
P = 0.320 days, and wc conclude our adopted 
spectroscopic period is not sensitive to the choice 
of the template. 

Adopting a period of P = 0.320 days, wc per- 
formed a finer search of the K2 — Tq piano, with 
step sizes of 1 km s^^ in K2 and 0.001 days in 
Tq. The results are shown in Fig. 4. The cross- 
correlation values have a fairly well-defined peak 
value along the Tq axis, where the maximum cross- 
correlation is for To = 2,452,436.600 (HJD). On 
the other hand, the peak along the K2 axis is 



broad with a maximum that occurs in the range 
430 < ^ 439 km s~^. Judging from the width 
of the peak and on our experience with HI 705-250 
(Remillard et al. 1996), we adopt K2 = 435 ± 30 
km s~^. The optical mass function is then 

PK^ 

/(M) = ^ = 2.73 ± 0.56 Mq. 

The best restframe spectrum found for the finer 
grid is shown in Fig. 5. Although few obvious K- 
star absorption features are apparent, the cross- 
correlation of this restframe spectrum using the 
K4V star BS 5568 as a template does yield a sig- 
nificant peak near zero velocity (Fig. 5). 

Since the spectroscopic period agrees with the 
photometric period, our initial determination that 
the K4V template provides the best match needs 
no refinement. In order to be more quantitative, 
we cross-correlated the best restframe spectrum 
from the fine grid against all six template spec- 
tra. Table 2 gives the cross-correlation value at 
the peak and the Tonry & Davis 'r' value (which 
is a measure of the significance of the peak) for 
each template. As noted earlier, the K4V tem- 
plate has the best cross-correlation, so we adopt 
a spectral type of K4V for the companion star. 
However, one should note that the GV5 star BS 
7330 provides the second-best match, and as such 
our assignment of a spectral type of K4V should 
be treated with caution. We attempted to decom- 
pose the restframe into its stellar and accretion 
disk components using the technique outlined in 
Marsh, Robinson, & Wood (1994). Using the K4V 
template BS 5568, we found that the accretion 
disk dominates in the P-band, where its contri- 
bution is > 80%. However, owing to the limited 
number of templates and the poor signal-to-noise, 
this result should also be treated with caution. 

Our measurement of the optical mass function 
sets a lower limit on the mass of the compact pri- 
mary. To find the actual mass we need to find the 
mass ratio of the two components and inclination 
of the orbital plane. If the companion star fills 
its Roche lobe and is in synchronous rotation, the 
mass ratio can be computed if one can measure K2 
and the projected rotational velocity of the com- 
panion star (e.g. Wade & Horne 1988). One can 
also get an estimate of the mass ratio if the radial 
velocities of the Ha emission line can be reliably 
measured (e.g. Orosz et al. 1994). However, in the 
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case of J1650, the poor quality of the spectroscopic 
data prevents us from making either measurement. 
We can use our Magellan light curve to place some 
constraints on the inclination. However, the re- 
sults are sensitive to the amount of light the ac- 
cretion disk contributes in the i?-band, and our 
estimate of the disk contamination outlined above 
is quite uncertain. Nevertheless, if we assume the 
limiting case of no light from the disk in the R- 
band, then we find an inclination of i = 50 ± 3° 
(see Fig. 3 for a representative model for this case). 
If we add a substantial contribution from the disk 
in the i?-band (> 80%), then the inclination is 
70° or higher (a representative ellipsoidal model 
for this case is also shown in Fig. 3). Owing to the 
lack of X-ray eclipses, the inclination is less than 
about 80° (the exact limit depends on the mass ra- 
tio). If, for the sake of discussion, we assume the 
mass ratio of J1650 is similar to that of A0620-00 
or GRS 1124-683 (Q w 10), then an inclination 
of 50° combined with our measured value of the 
optical mass function yields a black hole mass of 
7.3 ± 0.6 Mq. Likewise, an assumed inclination of 
70° yields a black hole mass of 4.0 ± 0.6 Mq. 

4. Discussion 

Is general relativity (GR) the correct theory 
of gravity in the strong fields found near a black 
hole? One promising approach to answering this 
fundamental question is offered by the key discov- 
ery of the Rossi X-ray Timing Explorer (RXTE) 
that seven Galactic black holes (including J1650) 
display quasiperiodic (QPO) X-ray oscillations in 
the range of 100 - 450 Hz (McClintock & Remil- 
lard 2004 and cited references). These fast QPOs 
must be produced near the event horizon since the 
X-rays originate there and since such high frequen- 
cies are comparable to the Kepler frequency of the 
innermost stable circular orbit around a black hole 
{vk = 2199 - 16150 (M/Mq)-i Hz depending on 
the dimensionless spin parameter a», where the 
full range of a* is to 1 e.g. Kato, Fukue, & Mi- 
neshige 1998). 

RXTE has made a further important discovery: 
four of these seven black holes produce pairs of sta- 
tionary, fast QPOs that have frequencies in a 3:2 
ratio (McClintock & Remillard 2004; Roman et al. 
2003a) . Such commensurate freqiiencies are a hall- 
mark of non-linear resonance phenomena. Thus, 



this discovery has promoted a "resonance model" 
that invokes enhanced emissivity at the radius in 
the accretion disk where two of the three spatial 
coordinate frequencies (e.g. Keplerian and radial) 
are commensurate (Abramowicz & Kluzniak 2001; 
Abramowicz et al. 2003). For the three black 
holes with the QPO pairs and measured masses 
(e.g. GRO J1655-40, XTE J1550-564, and GRS 
1915-1-105), the QPO frequency is correlated with 
the black hole mass, where v oc (e.g. McClin- 
tock & Remillard 2004). This scaling is expected 
for GR oscillations, but different sources can lie on 
the same curve only if they have similar values of 
the dimensionless spin parameter a* . 

Schnittman & Bertschinger (2004) performed 
ray tracing calculations in the Kerr metric for 
emitting blobs orbiting a black hole at a radius 
where there is a 3:1 resonance between the az- 
imuthal (fl) and radial coordinate frequencies. It 
was shown that general relativistic effects may im- 
part a QPO at and a beat oscillation at | f2 with 
a relative strength that depends on the angular 
width of the emitting blob of material. This in- 
terpretation, when combined with accurate mass 
measurements, yields values for the dimensionless 
spin parameter: 0.40 < a* < 0.55 for GRO J1655- 
40 and 0.3 < a* < 0.6 for XTE J1550-564 (Remil- 
lard et al. 2002). This illustrates the potential di- 
agnostic power of fast QPOs, if we can specify the 
correct oscillation mechanism. Furthermore, mea- 
surements of black hole mass and QPO frequency 
allow model-independent comparisons of spin dif- 
ferences between black hole subclasses as distin- 
guished, for example, by relativistic jets or binary 
period. 

With our measurement of the orbital period 
and optical mass function of J 1650, we have taken 
the first steps need to fully exploit the potcmtial 
of this system. If the QPO at = 250 ± 5 Hz 
observed by Homan et al. (2003b) represents the 
2i>Q oscillation, then the predicted mass of the 
black hole would be about 7.5 Mq, if J1650 lies 
on the same : M curve on which GRO J1655- 
40, XTE J1550-564, and GRS 1915+105 lie (i.e. 
if the black holes all have similar values of the 
dimensionless spin parameter a*). Although there 
are substantial uncertainties in both our measured 
value for the optical mass function and in our in- 
clination estimate, it appears that the black hole 
mass might be less than 7.5 Mq, which in the con- 
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text of the orbital resonance model would imply 
that the value a* in J1650 is much smaller than 
it is for the other three sources (e.g. Abramowicz 
et al. 2004). On the other hand, on the basis of 
X-ray spectroscopy, Miller et al. (2002) have ar- 
gued that J1650 contains a maximal Kerr black 
hole. Thus J1650 may be a challenging case for 
the orbital resonance model of the high frequency 
QPOs, although we again point out that our mass 
estimate has large uncertainties. 

Recent models for black hole formation (e.g. 
Fryer & Kalogera 2001) predict a continuous and 
roughly exponential distribution of masses for 
black holes in binary systems, where the num- 
ber of black holes formed falls off as the mass in- 
creases. Depending on their model assumptions, 
these authors predict the number of black holes 
with masses in the range 3 — 5 Mq should be 
roughly two to three times the number of black 
holes with masses in the range 8 — 10 Mq. Bailyn 
et al. (1998) performed a statistical analysis of 
seven black holes with dynamical mass estimates 
and found strong evidence for a peak in the mass 
distribution centered near 7 Mq and a "gap" in 
the distribution between the neutron star masses 
and the peak near 7 Mq. Since that time, the 
sample of objects has more than doubled, and it 
now appears that the evidence for a peak in the 
distribution near 7 Mq is weak (see the references 
cited in §1 and Orosz 2003). However, black holes 
with masses in the 3 — 5 Mq range still seem to 
be rare, with GRO J0422-I-32 being the only good 
candidate (Gelino & Harrison 2003). In this re- 
gard, J1650 is an interesting system because of its 
probable low mass. 

Clearly additional data will be required to make 
more definitive statements about the mass of the 
compact object in J1G50. It should not be un- 
duly difficult to obtain higher quality spectra. In 
hindsight, we now know that the companion star 
is a K-star, and as such it should have relatively 
strong absorption features near 5170A. This re- 
gion was not included in the VLT spectra owing 
to the relatively high resolving power of the 1200R 
grism. It should be easier to measure radial veloc- 
ities in spectra with more wavelength coverage (at 
the expense of resolving power). Also, spectral 
observations obtained in good to excellent seeing 
conditions (< 0.6 arcseconds) obviously will have 
better signal-to-noise than the current spectra do 



(the seeing varied between 1 and about 1.4 arcsec- 
onds). Additional photometry should be obtained, 
and it would be helpful if two or more bandpasses 
could be used, since the amount of contaminating 
light from the disk should vary with color. 

5. Summary 

Using i?-band photometry of XTE J1650-500 
obtained with the 6.5m Clay telescope, we have 
measured a photometric period of P = 0.3205 ± 
0.0007 days. This value is confirmed by our reanal- 
ysis of archival spectroscopic data. Our reanaly- 
sis of the spectroscopy also yields a velocity semi- 
amplitude for the companion star of K2 = 435 ±30 
km s~^, which when combined with the orbital 
period, gives an optical mass function of f{M) = 
2.73 ± 0.56 Mq. The spectral type of the compan- 
ion star is not well constrained because we only 
have six template spectra available to us. A tem- 
plate spectrum with a spectral type of K4V pro- 
vides the best match; next best matches are pro- 
vided by template spectra with spectral types of 
G5V and K2III, respectively. The summed spec- 
trum suggests that the accretion disk may con- 
tribute a substantial amount of light in R (« 80% 
or more). Using the -R-band light curve, we find a 
lower limit to the inclination of « = 50 ± 3° , which 
gives an upper limit to the mass of the black hole of 
Ml < 7.3 Mq, assuming a mass ratio of Q = 10. 
If we assume that the accretion disk contributes 
80% of the light in R, then the inclination would 
be about 70° or more, and the black hole mass 
would be about 4 Mq (again assuming a mass ra- 
tio of Q = 10). 

We thank Matt Holman and Paul Schechter 
for providing supplementary photometric data and 
Mauricio Navarette for help with the observations. 
This work was supported in part by NASA grant 
NAG5-9930. 
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Fig. 1. — An i?-band finding chart for J1650. Tlie field of view is one arcminute by one arcminute. 
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Fig. 2. — The results from our search for the binary period are presented. On the top, the periodogram 
made from the J1650 photometry using ELC is shown. This plot shows the absolute deviation as a function 
of the trial period, where a smaller value of the absolute deviation indicate a better fit. On the bottom, 
the periodogram for J1650 derived from the "restframe" analysis of the spectra is shown. This plot shows 
the maximum cross-correlation value as a function of the trial period, where a higher cross-correlation value 
indicates a better fit. Both periodograms have their best statistic at the same trial period, namely P = 0.320 
days (dot-dashed hne). The dotted line denotes the period reported by SF2002 {P = 0.212 days). 
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Fig. 3. — Top: The Magellan light curve of J1650 phased on the photometric period of P = 0.3205 days. 
Middle: The phased light curved smoothed using 20 bins in phase, along with two representative ellipsoidal 
model fits (see Table 1 for model parameters) The solid line shows a model with an accretion that dominates 
the light (denoted at model A in Table 1) and the dashed line shows a model with no accretion disk (model 
B in Table 1). Bottom: The Magellan light curved folded using the Sanchez-Fernandez period of P = 0.212 
days. 
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Fig. 4. — A grayscale image and contour plot showing the peak cross-correlation value as a function of K2 and 
To at a fixed period of P = 0.32 days. The maximum cross-correlation values happen when 430 ^ K2 ^ 439 
km s~^ and To = 2,452,436.600 (HJD). The thinner panels show the cross cuts of the cross-correlation 
values through the peak along the iir2-axis and the To-axis. 
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Fig. 5. — Bottom: The best normalized restframe spectrum of J1650 (top spectrum), the normalized spec- 
trum of the K4V star BS 5568, scaled by 0.1 and shifted upwards by 0.5 (bottom spectrum), and the 
difference spectrum (middle spectrum). Top: The cross-correlation function of the best restframe spectrum 
against the K4V star BS 5568. 
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Table 1 

Ellipsoidal Model Parameters 



Parameter'' Model A Model B 



i (deg) 


79.0 


51.7 


Q 


20.7 


1.1 


Ti (K) 


4500 


4500 




0.099 


0.099 


Idisk (K) 


15000 




dinner 


0.01 




^outer 


0.99 






-0.3 




Aim (deg) 


2.0 





•^See Orosz & Hauschildt (2000) for 
a full discussion of the parameters and 
their use in the ELC code. 



Table 2 

Cross-Correlation Values for Templates 



Star 


Spectral type 


Cross-Correlation 
Peak Value 


Tonry & Davis r 


BS 5568 


K4V 


0.36 


7.69 


BS 7330 


G5V 


0.27 


6.84 


HD 11301 


K2III 


0.26 


3.91 


BS 2668 


KOV 


0.18 


1.88 


BS 8042 


G3 IV 


0.10 


1.22 


HD 209290 


MOV 


0.06 


0.70 
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